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ABSTRACT 

We present the results of X-ray and near-infrared (NIR) observations of L1448, a star-forming 
region in the Perseus cloud complex using the Chandra X-ray Observatory and the 4 m telescope 
at the Kitt Peak National Observatory. We detect 72 X-ray sources in a '--^IT'xlT' region with a 
~68 ks ACIS exposure, for which we conduct follow-up NIR imaging observations in a concentric 
~ll'x 11' region with FLAMINGOS down to ~ 17 mag. Twelve X-ray sources have NIR or 
optical counterparts. By plotting X-ray mean energy versus NIR to X-ray flux ratio, the X-ray 
sources are clearly separated into two groups. The X-ray spectral and temporal features as well 
as NIR magnitudes and colors indicate that one group mainly consists of young stellar objects 
(YSOs) in the cloud and the other of background extragalactic sources. Ten X-ray-emitting YSO 
candidates are thus newly identified, which are low-mass or brown dwarf mass sources from their 
NIR magnitudes. In addition, a possible X-ray signal is found from a mid-infrared protostar 
L1448 IRS 3(A). The lack of detection of this source in our deep NIR images indicates that this 
source has a very steep spectral slope of >3.2 in 2-10 fim. 



Subject headings: X-rays: stars 
(L1448) 



infrared: stars — stars: pre-main-sequence — stars: individual 



1. INTRODUCTION 

L1448 is a dark cloud (Lynds 1962) located at 
the western edge of the Perseus molecular cloud 
complex at a distance of ^^300 pc. It has a 
size of ~15'x7' in a C^^O map and a mass of 
~100 Mq (Bachiller & Cernicharo 1986). The re- 
gion has been intensively studied in wavelengths 
longer than far-infrared focusing mainly on a gi- 
ant molecular outflow, and only few studies in the 
near-infrared (NIR) broad-band (Bally et al. 1993) 
and none in the X-ray band have been presented, 
both of which are currently strong measures to 
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identify stellar constituents of star-forming clouds. 
The infrared excess due to emission from circum- 
stellar material and the X-ray brightness are two 
major observational characteristics of young stel- 
lar objects (YSOs), which can be respectively re- 
vealed through NIR broad-band photometry in the 
J, H, and Ks bands and X-ray imaging observa- 
tions. The purpose of this paper is to report the re- 
sults of the first X-ray imaging-spectroscopy study 
of LI 448 aided by follow-up NIR broad-band imag- 
ing observations. Previous X-ray studies in the 
Perseus cloud complex can be found in Preibisch 
et al. (1996); Preibisch & Zinnecker (2001, 2002, 
2004) for IC348, Getman et al. (2002); Preibisch 
(2003) for NGC 1333, and Yamauchi et al. (2001) 
for Barnard 1 and IRAS 03282-1-3035. 

The most magnificent feature in LI 448 is a 
parsec-scale bipolar outflow emanating from the 
core of the cloud, which is one of the highest- 
velocity, best-coUimated, and youngest outflows 
known (Davis & Smith 1995). The outflow was 
discovered by Bachiller et al. (1990) in a molec- 
ular emission line image. The exciting source of 
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the outflow was promptly detected in the millime- 
ter and centimeter continua (Bachiller et al. 1991; 
Curiel et al. 1990), which is referred as L1448 
mm or L1448 C(entcr). Barsony et al. (1998) ob- 
tained a spectral energy distribution (SED) in far- 
infrared to millimeter wavelengths and established 
the class status of this source. 

The Infrared Astronomical Satellite (IRAS) de- 
tected three sources (IRAS 03220-^3035, 03222+3034, 
and 03225+3034) in this region, which are called 
L1448 IRS 1, 2, and 3. IRS 1 has a visual counter- 
part with mv = 19.5 mag (Cohen & Kuhi 1979). 
From its strong Ha emission and excess emission 
in the IRAS band (Cohen & Kuhi 1979; de Grijp 
et al. 1987), it is probably a classical T Tauri star 
or a class I protostar in the cloud. IRS 2 is a 
class source confirmed by O'Linger et al. (1999) 
based on an SED using IRAS and SCUBA data. 
This source is associated with various signatures 
of extreme youth, including CO and II2 molecu- 
lar outflows (O'Linger et al. 1999; Eisloffel 2000), 
Herbig-Haro objects (Bally et al. 1997), an H2O 
maser source (Anglada et al. 1989), and a centime- 
ter continuum emission (Anglada & Rodriguez 
2002). IRS 3, also known as L1448 N(orth), is 
resolved into three components; IRS 3(A), 3(B), 
and 3(C), also known respectively as L1448 N(A), 
N(B), and NW. The former two constitute a bi- 
nary separated by ~ 2100 AU. Classification of 
the individual binary components in millimeter 
to far-infrared observations had been difficult due 
to crowding. However, recent mid-infrared (MIR) 
observations clearly resolved IRS 3(A) and 3(B) 
with an iV-band detection only from 3(A). Based 
on this, Ciardi et al. (2003) argued that 3(B) 
is a class protostar whereas 3(A) might be at 
the transition phase between class and I stages. 
IRS 3(C), away from the binary by ~20", is an- 
other class source (Barsony et al. 1998). 

2. OBSERVATIONS 

The X-ray observation was carried out on 
2004 March 15 using ACIS (Advanced CCD 
Imaging Spectrometer; Garmire et al. 2003) on- 
board Chandra (Weisskopf et al. 2002). Four 
ACIS-I chips (10, II, 12, and 13) covered a 
-17'xl7' region aimed at R. A. = 3^25™36.4« and 
Decl. = 30°44'58" (.12000.0) for an exposure time 
of 67.9 ks. ACIS has sensitivity in the 0.5-9.0 keV 



band with a resolution of ^150 eV at 6 kcV and 
a superb point spread function (PSF) radius of 
~0.5" at the on-axis position. The data were 
taken with the very faint telemetry mode and the 
timed exposure CCD operation with a frame time 
of 3.2 s. 

The NIR observation was conducted on 2003 
December 16 using FLAMINGOS (Elston et al. 
2003) on the Cassegrain focus of the 4 m telescope 
at the Kitt Peak National Observatory. FLAMIN- 
GOS uses a HgCdTe HAWAII-2 array which has 
a format of 2048x2048 with a pixel scale of 
~0.316" pixel^^. With dithering observations of 
an amplitude of 15", we obtained J-, H-, and Ks- 
band images that cover a ~ll'xll' area centered 
close to the ACIS aim point at R. A. = 3'^25'"35.1" 
and Decl.=30°45'19" (J2000.0). The total expo- 
sure time was 13.5, 7, and 17 minutes respectively 
in the J, H, and Ks bands. The seeing was ^1". 
Figures 1 (a) and (6) show the X-ray and NIR 
intensity maps. The outflow from L1448 mm ap- 
pears clearly in the Ks band (Fig. 16) presumably 
due to II2 v=l-0 S(l), vibrational-rotational tran- 
sition line emission by excited hydrogen molecules. 

3. DATA REDUCTION & ANALYSIS 
3.1. X-ray Data 

We reprocessed the Level 1 data distributed by 
the Chandra X-ray Center to obtain our X-ray 
event list. The latest calibration results (CALDB 
2.28) were incorporated and the backgroimd re- 
jection algorithm specific to data taken with the 
very faint mode was applied. Events were further 
cleaned by removing cosmic-ray afterglows and ap- 
plying filters based on the event grade, status, and 
good-time intervals. 

We detected sources using a wavelet technique 
in CIAO^ from three images of different energy 
bands; soft (0.5-2.0 keV), hard (2.0-9.0 keV), and 
total (0.5-9.0 keV). For each detected source, we 
used ACIS Extract^ for a systematic event extrac- 
tion, background subtraction, calculation of in- 
strumental responses, and construction and bin- 
ning of spectra and light curves. Source events 
were extracted from a 90% encircled energy poly- 

^Scc http://asc.harvard.edu/ciao/ for detail. 

^Scc http://www.astro.psu.edu/xray/docs/TARA/afi_tisers_guide.html 
for detail. 
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Fig. 1. (a) X-ray image. The 0.5-9.0 keV in- 
tensity is given by a logarithmic grey scale. The 
positions of X-ray sources with and without NIR 
counterparts are shown by crosses and pluses re- 
spectively, whereas the position of three IRAS 
sources (IRS 1-3) and L1448 mm are by the filled 
squares. The contours indicate the intensity of an 
NHs transition line (Anglada et al. 1989). The 
fields of view of four ACIS-I chips and FLAMIN- 
GOS are illustrated with the solid oblique squares 
and the dotted square. (6) NIR image. The Ks- 
band intensity is given by a logarithmic grey scale. 
The positions of embedded sources are bracketed 
with a pair of solid lines when undetected and by 
crosses when detected. 



gon of the 1.5 keV PSF, whereas background 
events were locally extracted from an annular re- 
gion around the source. The positions of sources 
were fine-tuned by correlating the event distribu- 
tion with the PSF. As the X-ray source distribu- 
tion is relatively sparse in the study field (Fig. la), 
we do not suffer any overlaps of extraction regions. 
The details of the procedure are described in Get- 
man et al. (2005). 

For each source, wc; derived X-ray photometry 
information consisting of the source position and 
its uncertainty, source count (0.5-8.0 keV), net 
comit rate (NCR), probability of no source (PNS), 
photometry significance (PS), and mean energy 
(ME). The net count (Cnet) is calculated from the 
number of counts (Csrc and Cbkg) and the extrac- 
tion area (Agrc and ^bkg) of the source and back- 
ground regions as Cnct — C*src — ^bkg^ (^src/^bkg), 

which is divided by the effective exposure time to 
derive NCR. PNS is an index for the detection sig- 
nificance, with which we test the null hypothesis 
that all detected counts are explained by back- 
ground fluctuation; 



PNS 
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diVP(Cbkg 
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bkg 



.^), (1) 



where P(A, N) represents a Poisson probability 
distribution function of a mean A to have N 
counts. PS is a metric for the reliability of pho- 
tometry defined as 

pg ^net 

~ v'(AC,rc)2 + (ACbkg)2 X (^rcMbkg)2 ' 

(2) 

where the approximation by Gehrels (1986) is 
used to estimate ACsrc 

1 + 



1 + y/C'src + 0.75 and 
ME is a metric to 



ACbkg = 1 + yCbkg + 0.75 
represent the spectral hardness of sources. The 
value is the average incident energy (E) weighted 
by Cnet of each energy bin and is defined as 



ME: 



jdEEC^,,{E) 
JdECnetiE) ■ 



(3) 



Table 1 shows the results of the ACIS photome- 
try, consisting of 72 sources with high significance 
both in the detection (PNS < 1x10"^) and the 
photometry (PS > 2). 

For each source, flux variability was tested us- 
ing the K-S test and the Bayesian block segmenta- 
tion technique (Scargle 1998). We recognized flux 
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Table 1 
ACIS Photometry 



x-ray* 


R.A.*" 


Decl.'" 


Error 


Counts 


NCR" 


PNS° 


PS" 


ME" 


NIR"* 


ID 


(hh:mm:ss.s) 


(dd:inm:ss) 


(") 




(ks-i) 






(keV) 


Counterpart 


1 


3:24:59.3 


30:44:41 


0.7 


25 


0.3 





3.2 


3.6 




2 


3:25:02.7 


30:38:34 


1.1 


25 


0.6 


9c-7 


2.8 


3.3 




3 


3:25:03.8 


30:43:28 


0.6 


20 


0.3 





3.0 


2.7 




4t 


3:25:04.8 


30:38:58 


0.4 


116 


2.4 





9.0 


1.4 


03250474+3038578 


5 


3:25:07.2 


30:38:04 


1.0 


34 


0.6 





3.8 


3.6 





^Variable sources are marked with a dagger. 
''Equinox in ,12000. 0. 

"See text for definitions of NCR, PNS, PS, and ME. PNS values are replaces with when PNS < le-9. 

*'2MASS counterparts are given in the format of hhmmssss+ddmmsss, while FLAMINGOS counterparts are given 
by their source numbers in Table 3. 

Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains 
only a sample. 



to be variable if the number of Bayesian blocks ex- 
ceeds two and the K-S probability is less than 1%. 

Four sources (Nos. 4, 7, 34, and 69) have variable 
flux, which are marked with daggers in Table 1. 

For 18 sources with more than 40 counts, we 
performed spectral analysis. The spectra were 
fit with a one-temperature thin-thermal plasma 
model (the APEC model; Smith et al. 2001) with 
the metallicity value fixed to 0.3 solar to deter- 
mine the best- fit absorption (A^h), plasma tem- 
perature (fceT'), and X-ray luminosity (Lx) val- 
ues. They were also fit with a power-law model to 
derive the best-fit Ah, index (F), and absorption- 
corrected flux {F^'^) values. Table 2 shows the 
results of 13 and 11 successful fits with the ther- 
mal and power-law models with a reasonable value 
of kBT< 10 keV and F <3. 

3.2. NIR Data 

All FLAMINGOS frames were reduced follow- 
ing the standard procedures using IRAF; i.e., sub- 
traction of dark current, flat-fielding using a dome 
flat, subtraction of sky using median sky, removal 
of bad pixels, and trimming of edges. Several 
frames with readout failure were discarded. Arti- 
ficial signals were masked out, which we found at 
pixels to the right and left of, or above and below 
bright sources by ^128 pixels. Nebulosity extend- 
ing larger than the dithering amplitude may re- 
main in the median-sky-subtracted images; it does 
not affect our photometry. Referring to sources in 



the 2MASS all-sky survey catalog (Skrutskie et al. 
1997), an FLAMINGOS frames were corrected for 
astrometry by a low-order polynomial function to 
compensate for the image distortion. They were 
also corrected for photometry to match 2MASS 
magnitudes. After these corrections, the frames 
were combined into the final J-, H-, Ks-hand im- 
ages. 

We extracted sources from the Ks-hand im- 
age (Fig. 16) above a 5 a level using SExtractor 
(Bertin fc Arnouts 1996). The number of false 
positive det(X'tions is statistically negligible. Af- 
ter manually removing suspicious sources (intro- 
duced mostly by misidentifying outflow knots as 
point sources) and sources at the edges, we ob- 
tained 294 detections. We measured J-, H-, and 
Ks-ha.nd magnitudes at their positions (Table 3). 
We matched the FLAMINGOS and the 2MASS 
sources and found that all the 67 2MASS sources 
in the FLAMINGOS fleld of view have FLAMIN- 
GOS counterparts. 

Figure 2 shows the astrometric accuracy of 
the FLAMINGOS sources, where the difference 
of FLAMINGOS and 2MASS positions (AR. A. 
and ADecl.) are plotted for the 67 FLAMINGOS- 
2MASS counterpart pairs. The root mean square 
of AR. A. and ADccl. are 0.18" and 0.33", indi- 
cating that the FLAMINGOS position accuracy 
is ~0.4" («1.2 FLAMINGOS pixel). The accu- 
racy is better for most sources, as the values are 
exaggerated by three sources at the top left in Fig- 
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Table 2 

ACIS Spectroscopy of Bright Sources 







— thermal model 


















power-law model 


X-ray 








JVh" 


pa 




ID 


(10^^ cm-2) 


(keV) 


(10" ergs s-^) 


(10^^ cm- = ) 




(10-^* ergs s-^ cm-^) 


4 


0.0 (0.0-0.1) 


1.7 (1.5-1.9) 


1.8 (1.1-2.0) 








6 


0.4 (0.2-0.6) 


0.9 (0.7-1.1) 


1.3 (0.0-1.6) 








10 


1.3 (0.8-2.0) 


2.7 (1.4-5.9) 


3.6 (0.2-3.8) 


1.6 (1.3-2.3) 


2.7 (2.3-3.3) 


5.7 (0.0-6.7) 


11 








0.3 (0.2-0.6) 


1.2 (0.9-1.5) 


4.1 (0.0-4.6) 


15 


0.0 (0.0-0.2) 


1.0 (0.9-1.2) 


1.0 (0.0-1.5) 








16 


0.3 (0.2-0.4) 


5.3 (4.3-7.1) 


24.9 (21.7-26.5) 


0.4 (0.3-0.5) 


1.9 (1.8-2.0) 


26.1 (21.7-29.0) 


20 


0.3 (0.2-0.3) 


7.9 (5.5-11.2) 


10.2 (8.3-10.7) 


0.3 (0.3-0.4) 


1.7 (1.6-1.8) 


10.1 (7.8-11.4) 


23 


0.3 (0.1-0.5) 


0.6 (0.5-0.8) 


0.9 (0.4-1.3) 








42 


0.0 (0.0-0.1) 


0.6 (0.5-0.6) 


2.1 (0.8-2.0) 








48 


0.9 (0.7-1.0) 


0.6 (0.5-0.6) 


7.7 (3.6-12.4) 








50 








0.1 (0.0-0.3) 


1.4 (1.0-1.6) 


1.3 (0.0-1.5) 


53 








0.4 (0.1-0.9) 


1.1 (0.6-1.6) 


1.4 (0.0-1.6) 


54 








0.5 (0.3-1.0) 


1.6 (1.2-2.1) 


1.7 (0.0-2.1) 


64 


0.8 (0.6-1.0) 


5.7 (3.8-9.9) 


7.1 (3.8-7.4) 


0.9 (0.7-1.2) 


1.9 (1.7-2.1) 


7.5 (3.6-8.8) 


66 


0.4 (0.2-0.7) 


3.2 (2.0-6.4) 


1.3 (0.3-1.7) 


0.5 (0.3-0.9) 


2.3 (1.8-2.9) 


1.5 (0.0-1.9) 


68 


2.1 (1.4-3.1) 


2.0 (1.2-4.7) 


3.9 (0.2-4.3) 








71 


0.8 (0.3-1.4) 


5.7 (2.3-15.0) 


3.3 (0.0-3.7) 


0.8 (0.5-1.6) 


1.8 (1.3-2.7) 


3.3 (0.0-4.3) 


72 








1.0 (0.8-1.4) 


1.6 (1.3-1.9) 


5.1 (0.0-5.9) 



^The 1 fj confidence ranges are given in the parentheses. 

^Absorption-corrected values in the 0.5—8.0 keV band. A distance of 300 pc is assumed to derive Lx for the thermal fits. 



Table 3 
FLAMINGOS Photometry 



NIR" 


R.A.'' 


Docl.'' 


./'= 




Ks 


2MASS 


ID 


(hh:mm:ss.s) 


(dd:mm:ss) 


(mag) 


(mag) 


(mag) 


Counterpart 


1 


3:25:09.4 


30:48:59 


18.8 


18.1 


17.1 




2 


3:25:09.5 


30:46:21 


12.5 


11.0 


10.1 


03250943-1-3046215 


3 


3:25:09.7 


30:48:48 


13.3 


12.5 


12.2 


03250965-1-3048488 


4 


3:25:10.1 


30:46:07 


16.4 


16.2 


17.6 




5 


3:25:10.2 


30:47:13 






16.5 





^Sources with NIR excess are marked with a dagger. 

^Equinox in J2000.0. 

'^Magnitudes for sources detected with more than 3 a are given. 



Note. — The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 
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urc 2, all of which are located too far away from 
the image center to be corrected for distortion. 



o 

(D 

O 



o 

Q 
<1 



0.5 




0.0 



-0.5 - 



-1.0 -0.5 0.0 0.5 
AR.A. (arcsec) 



1.0 



X-ray sources were identified in the optic;al c;ata- 
log. No X-ray emission was found from Herbig- 
Haro objects (Bally et al. 1997) or H2 knots (Davis 
& Smith 1995). 

4. DISCUSSION 

4.1. Nature of X-ray Sources 

Two independent quantities from the X-ray 
photometry (ME and NCR) combined with the 
NIR identification give a clue to classify and re- 
veal the nature of the 72 X-ray sources. Figure 3 
shows a scatter plot of ME and NCR. Most of the 
X-ray sources with NIR counterparts {filled sym- 
bols) are distributed in the limited ME range of 
<2.5 keV, while those without NIR counterparts 
{open symbols) are distributed in the limited ME 
and NCR ranges of >2.0 keV and <1.0 ks'^ This 
implies that the X-ray sources consist of two major 
classes of different nature. 



Fig. 2. — Astrometric accuracy of FLAMINGOS 
sources. The difference of R. A. and Decl. be- 
tween the FLAMINGOS and 2MASS positions 
(AR.A. and ADecl.) are plotted for all the 
FLAMING0S-2MASS pairs. The solid square at 
the center represents the FLAMINGOS pixel size 
of 0.316"x0.316". 

We estimated the detection completeness of 
the FLAMINGOS images by embedding artificial 
sources and detecting them with the same source 
detection algorithm. The detection rates at differ- 
ent magnitude bins were derived to estimate the 
90% completeness limit as ^18.0 mag, ~^16.5 mag, 
and ^^17.0 mag in the J, H, and Kg bands. 

3.3. NIR and Optical Identification of X- 
ray Sources 

Seven X-ray sources were found to have FLAMIN- 
GOS counterparts within 1.5" and five additional 
sources to have 2MASS counterparts outside of the 
FLAMINGOS view. X-ray positions in Table 1 are 
shifted by ~0.3" to match the 2MASS astrometry. 
The root mean square of the positional difference 
between 2MASS-X-ray pairs is -0.6" («1.3 ACIS 
pixel). Ten of the NIR- identified X-ray sources 
are optically identified with the USNO-B catalog 
(Monet et al. 2003). None of the NIR- unidentified 
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o 




12 3 4 5 
Median Energy (keV) 



Fig. 3.— Scatter plot of ME and NCR for all X- 
ray sources. Sources with and without NIR coun- 
terparts are shown with filled and open symbols, 
respectively. The source numbers (Table 1) of the 
former are labeled in Roman, while the latter is 
in Italic if they have >40 counts. The median 
value of ME and NCR are shown with dashed- 
and-dotted lines. 

The separation of the two classes becomes more 
apparent when NIR flux is taken into account. 



6 



Figure 4 shows a scatter plot of ME and the ra- 
tio of NIR fliix in the Kg band to X-ray flux in 
0.5-8.0 keV band. Here, the X-ray flux (Fx) is 
calculated by 



Fx = 



JdE E Cr^etiE) _ ME X NCR 



EA t. 



exp 



EA 



(4) 



where te^p is the exposure time and EA is the 
effective area of the mirror averaged over the en- 
ergy. The X-ray sources with NIR counterparts 
{large filled symbols) and those without FLAMIN- 
GOS counterparts {large open squares) are clearly 
separated into two groups. One group (group A) 
is comprised of X-ray source Nos. 4, 6, 7, 14, 15, 
23, 30, 42, and 48 in the upper left quarter, while 
the other (group B) includes Nos. 11, 16, 39, and 
FLAMINGOS-unidentified X-ray sources {large 
open squares) at the bottom. NIR- unidentified 
X-ray sources outside of the FLAMINGOS field 
of view {small open symbols) are hard to classify 
into either group as the Ks-hand upper limit of 
2MASS is not deep enough. This illustrates the 
importance of our NIR observation deeper than 
2MASS, even if it results in no detection from 
X-ray sources. 

We conclude that the sources in group A arc 
mostly YSOs that belong to L1448, while those in 
group B are background AGNs. We note that NIR 
identification of an X-ray source in star-forming 
regions does not always guarantee its stellar ori- 
gin because AGNs, which constitutes the largest 
contaminant, can be detected even at the 2MASS 
depth. Instead, we employ the NIR to X-ray flux 
ratio as an effective measure to separate these two 
classes. The optical to soft X-ray flux ratio was 
used for a similar classiflcation purpose in X-ray 
blank-sky surveys by Einstein and ROSAT satel- 
lites (Maccacaro et al. 1988; Krauttcr ct al. 1999). 
We use the NIR to total band X-ray flux to be suit- 
able for star-forming regions, which are routinely 
highly obscured. In (low-mass) main sequence and 
pre-main-sequence (PMS) sources. X-ray emission 
originates from stellar coronae while NIR emission 
originates from photosphere and from circumstcl- 
lar matter. In AGNs, on the other hand. X-ray 
emission is from accretion disks around the cen- 
tral black hole while NIR is from hot dust and jets. 
The different origins bring different typical values 
of NIR to X-ray flux ratio in stars and AGNs. 

Several lines of evidence support this conclu- 
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Fig. 4. — Scatter plot of ME and the flux ratio 
of NIR {Ks) and X-ray (0.5-8.0 kcV) bands for 
all X-ray sources with and without NIR counter- 
parts {filled and open). Open symbols have cither 
a normal size (for sources inside the FLAMINGOS 
field of view) or a smaller size (otherwise). Vari- 
able sources arc marked with diamonds. Open 
symbols indicate the upper limit arising from no 
Kg-hand detection. Labels are given in the same 
manner with Fig. 3. Variable sources are also la- 
beled. The median values of ME is shown with 
the dashed-and-dotted line. 
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sion. First, all sources with >40 counts in group 
A (Nos. 4, 6, 15, 23, 42, and 48) have a success- 
ful fit with thermal spectral models (Table 2) but 
not with power-law models with reasonable pa- 
rameters. On the other hand, all sources with 
>40 counts in group B (Nos. 11, 16, 50, and 53) 
have a successful fit with power- law models. Three 
(Nos. 11, 50, and 53) of them are not fit with 
thermal models (Table 2). The clear distinction 
in X-ray spectral features is consistent with the 
idea that (1) these groups of sources are different 
in nature, (2) group A sources are stellar X-ray 
emission, which typically show thermal spectra of 
ksT^ 1 keV, and (3) group B sources are AGNs, 
which typically show power-law spectra of F ~1.7. 
This conversely enables us to use X-ray spectral 
features to classify sources for which it is unclear 
which group to belong to. Nos. 54 and 72 are group 
B sources for successful fits only by the power-law 
model (Table 2). 

The second argument is the X-ray flux variabil- 
ity of a typical time scale of ~10 ks. Stellar X-ray 
emission from coronae is often characterized by X- 
ray flares due to magnetic reconnection, whereas 
no such short time-scale variability is expected 
from AGNs. Sources with X-ray variability are 
seen in group A (Nos. 4 and 7) and not in group B. 
This conversely indicates that sources with similar 
variability (Nos. 34 and 69) belong to group A. 

Third, NIR magnitudes of the X-ray sources in 
group A are distributed in ttiks ^ 9-13 mag. Con- 
tamination by extragalactic sources in this magni- 
tude range is negligible as the expected number of 
galaxies in the FLAMINGOS field of view is <1 
in m.Ks ^ 14 mag (Tokunaga 2000). The color- 
magnitude diagram (Fig. 5) shows that group A 
sources have consistent NIR magnitudes for low 
mass or very low mass stars at a distance of 300 pc. 
Sources in group B, on the contrary, have uiks ^ 
15 mag (Tables 1 and 3), which is too faint to be 
stars at that distance. 

Fourth, the NIR to X-ray flux ratio of group 
A sources {Fk^/Fx ^ 10^-10^) is consistent with 
typical values of the X-ray to the bolometric lumi- 
nosity ratio of Lx/L^oi ~ 10~^-10~^ (Preibisch 
et al. 2005) of PMS sources. Assuming an age 
of 2 Myr (the same with IC 348; Muench ct al. 
2003), and using Lxs/Lboi ~ 4-9x10-^ (0.08-1.4 
Mq PMS source; Baraffe et al. 1998; D'Antona 
& MazzitelU 1994) and Lks/Lx ~ Fks/Fx ~ 
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Fig. 5. — Color-magnitude diagram of X-ray 
sources with NIR counterparts {filled squares) and 
FLAMINGOS sources {open circles). Sources 
with J-, H-, and Ks-hand detections with magni- 
tude uncertainty of < 0.05 mag are plotted. Labels 
are given to the X-ray sources with NIR counter- 
parts. The 1, 2, and 3 Myr isochrone curves {dot- 
ted, solid, and dashed curves) are from Baraffe 
et al. (1998). The reddening lines for 1.4 Mq, 
0.2 Mq, and 0.08 Mq are given with dashed lines 
assuming an age of 2 Myr. 
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5x10^ (the median value of group A sources), we 
obtain the value of Lx/^boi ~ 0.8-1.8x10"'* for 
group A sources. A similar Lyi/LhoX value is ob- 
tained if we assume that these sources arc main se- 
quence stars (Tokunaga 2000; Drilling & Landolt 
2000), which is higher than their typical values 
(IQ-'^ 10-'^ Preibisch et al. 2005). We thus con- 
sider that most of the group A sources are YSOs 
in the cloud and not main sequence stars in the 
line of sight. The distinction between YSO and 
main sequence natures should be confirmed with 
follow-up spectroscopic studies. Group B sources, 
on the other hand, have Fks/Fx values of 10~* 
10^, which is consistent with the values seen in an 
X-ray-selected sample of AGNs (Watanabe et al. 
2004). 

Figure 6 shows a color-color diagram of NIR- 
identified X-ray {squares) and FLAMINGOS (cir- 
cles) sources. All NIR-identified X-ray sources in 
group A (Nos.4, 6, 7, 14, 15, 23, 30, 42, and 48) 
have NIR colors consistent with being reddened 
weak-line T Tauri or main sequence stars. None of 
them shows apparent NIR excess emission, which 
does not exclude the possibility that they are clas- 
sical T Tauri stars or protostars, as Ks-hand ob- 
servations are not completely sensitive to excess 
emission from circumstellar disks (e.g., Haisch et 
al. 2001). The only source with significant J-, H-, 
and Ks-ha,nd detections in group B (No. 16) has an 
inconsistently red color for reddened stars of mj— 
rriKs = 2.0 mag, which is commonly seen among 
2MASS red AGNs (Cutri et al. 2001). The density 
of such sources is ~0.5 arcdeg"^ at Ks < 14.5 mag 
(Cutri et al. 2001), making the probability ^5% 
to have a 2MASS red AGN in an ACIS field. 

In total, twelve X-ray sources (Nos. 4, 6, 7, 14, 
15, 23, 30, 34, 42, 48, 68, and 69) have X-ray and 
NIR features consistent of being low-mass YSOs in 
this cloud. Two of them are known sources; No. 7 
is IRS 1 and No. 42 is an A5 type star (Roeser & 
Bastian 1988) that consists of a double (Worley & 
Douglass 1997). Main sequence A type stars are 
not expected to have X-ray emission due to the 
lack of any plasma creation mechanisms (Schmitt 
et al. 1993). The X-ray emission in No. 42 may 
come from a low-mass companion, which mim- 
ics a I Fx value consistent with a PMS source 
(Fig. 5). The remaining ten sources are candidates 
of new YSOs in L1448. The brightness histogram 
of the rest of the X-ray sources is consistent with 



a log TV-log S relation of AGNs (Moretti et al. 
2003). 
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Fig. 6. — Color-color diagram of X-ray sources 
with NIR counterparts (filled squares) and 
FLAMINGOS sources (open circles). Sources 
with J-, H-, and i^g-band detections with mag- 
nitude uncertainty of < 0.1 mag are plotted. 
FLAMINGOS sources with NIR excess emission 
are shown with larger symbols. The reddening 
lines are shown (dashed lines) from the intrinsic 
colors of dwarfs and giants (fMck solid curves; 
Tokunaga 2000) and the classical T Tauri star lo- 
cus (thick solid lines; Meyer et al. 1997). 

4.2. X-ray and NIR Emission from Em- 
bedded Sources 

4-2.1. Infrared Excess Sources 

From Figure 6, ten FLAMINGOS sources (large 
open circles) are found to have NIR excess emis- 
sion and thus are candidates for new embedded 
sources. Daggers are given for these sources in 
Table 3. Among them, FLAMINGOS Nos. 23 and 
149 are promising YSO candidates with estimated 
magnitude uncertainty of < 0.05 mag. 

4.2.2. IRS 1 

Both X-ray and NIR emission was found from 
IRS 1 (ACIS No. 7 in Table 1 and FLAMINGOS 
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No. 2 in Tabic 3 or 2MASS 03250943+3046215). 
The X-ray and NIR properties are consistent with 
the nature of this source as a classical T Tauri 
star or a class I protostar; (1) it shows X-ray 
flux variation commonly seen in X-ray emission 
of stellar magnetic activity origin, (2) its ME and 
Fks/Fx values indicate this source to be a YSO 
(Fig. 4), and (3) a nebulous structure is seen in the 
FLAMINGOS Ks-hand image (Fig. 16). The NIR 
magnitudes of this source (toj ^ 12.5 mag, ttih ^ 
11.0 mag, and m^s ~ 10.1 mag) correspond to 
an object with the mass of 0.5-1.0 Mq at 300 pc 
obscured by Ay 7 mag. 

4.2.3. IRS 3(A) 

No emission is seen at the position of IRS 3(A) 

in our NIR images. The 5a level of the local sky 
noise corresponds to ~18.8 mag, ~17.7 mag, and 
~17.5 mag in our J-, H-, and Ks-h&nd images. 
The SED slope defined as a = (d log XFx ) / (d log A) 
between 2.2 /zm {Ks band) and 10 /xm {N band) 
is commonly used to classify YSOs (Wilking et 
al. 1989), using the fact that the 2-10 ^m SEDs 
at a younger evolutionally stage are more dom- 
inated by emission from circumstellar material. 
The lower limit of a of IRS 3(A) is >3.2 from the 
non-detection in the NIR bands and a MIR de- 
tection in the N band (Ciardi et al. 2003). Here, 
we calculated the a value as the index between 
H and A'^ bands, as the ifg-band image is noisy 
at the position due to the emission arising from 
the outflow of IRS 3(A) itself. The value thus is 
a lower limit also in the sense that if-band emis- 
sion is subject to contamination from reflected 
photospheric emission. 

Greene et al. (1994) gave a qualitative classi- 
fication of YSOs using a, where a < —1.6 for 
class III, —1.6 < a < —0.3 for class II, and 
a > —0.3 for class I sources. Almost all known 
class I protostars have a < 2 with the only excep- 
tion of WL 22 {a = 3) in /9 Ophiuchi dark cloud 
(Wilking et al. 1989). Sources with a steeper in- 
dex than ~ 3 have been rarely known, as they are 
quite difficult to detect in K band images that tra- 
ditionally pilot protostar searches. However, re- 
cent progress of sensitive MIR imaging both by 
ground- and space-based telescopes is changing 
the situation by reporting several protostars with 
very steep spectra with indices of a > 3 (Cep 
E by Noriega-Crespo et al. 2004 and source Xe 



in R CiA by Hamaguchi et al. 2005 and private 
communications with K. Nedachi). The index of 
>3.2 of IRS 3(A) places this source in the class of 
very steep spectrum protostars. More sources in 
this emerging class are expected to be discovered 
with Spitzer observations of near-by star-forming 
regions. Some of them are class protostars (e.g., 
Cep E mm; Noriega-Crespo et al. 2004) and oth- 
ers are iiT-band- unidentified class I protostars, or 
in between. 

An unusual concentration of three X-ray pho- 
tons within a ~0.5" radius (~90% encircled en- 
ergy radius at -^0.4' off-axis angle) can be found 
at the position of IRS 3(A). No such concentration 
is seen within 1' of the on-axis position. Although 
it was not recognized as a source in the initial X- 
ray source search using the wavelet technique, its 
PNS value of 5x 10~* indicates the existence of an 
X-ray source of Fx ~ 5x10""'^^ ergs cm~^. The 
facts that the position of the concentration is con- 
sistent with that of IRS 3(A), that the energy of 
all photons is concentrated in the hard (2 5 keV) 
band, and that the arrival times of all photons fall 
in a limited time range of ~15 ks in the ~68 ks 
observation, which is indicative of a flare, may fur- 
ther strengthen the claim. There is no indication 
of temporal non-uniformity of background counts. 
A similar weak but plausible X-ray emission was 
found in an embedded source (FIR4) in NGC 2024 
by Chandra (Skinner et al. 2003). 

Given the capability of Chandra and MIR tele- 
scopes (Spitzer and ground-based facilities), it is 
now becoming possible to investigate X-ray emis- 
sion from protostars with a very steep spectrum 
(a > 3) even without Ks-hand identifications. If 
IRS 3(A) is indeed an X-ray source, it has the 
steepest spectral index among YSOs with X-ray 
emission. Although it is not yet certain whether a 
steeper a value indicates a younger age of pro- 
tostars, X-ray observations for steep spectrum 
sources are a key to answer how early X-ray emis- 
sion starts and how intense it is at the earliest 
stage of star formation. These are important ques- 
tions to understand star formation because the 
coupling between gas and magnetic fields can be 
largely afi^ected by X-ray ionization (Feigelson & 
Montmerle 1999). 
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4-2.4- Class Sources 

No significant X-ray or NIR emission was found 
from the four known class sources; L1448 mm, 
IRS 2, IRS 3(B), and IRS 3(C). The upper limit 
(PNS = 1 X 10"'^) of the X-ray flux in the 0.5 
8.0 keV is ~0. 5-1.0x10"^^ ergs s cm~^, whereas 
the upper limits (5 a) in the NIR regime are 
~18.8 mag (J), ^17.7 mag {H), and ^17.5 mag 
{Ks), except for L1448 mm and IRS 3(B) in the 
Ks-ha.nd, which are contaminated by local diffuse 
emission. 

The lack of X-ray detection from class sources 
is widely reported. Froebrich (2005) compiled 28 
well documented bona-fide class sources. Six- 
teen of them were observed by Chandra and nine 
in references are reported to have no X-ray detec- 
tion (L1448 NW and C; this work, NGC 1333 12, 
I4A, and I4B; Getman et al. 2002, OMC 3 MMS6; 
Tsujimoto et al. 2002, HH24 MMS and HH25 
MMS; Simon et al. 2004, and VLA 1623; Imanishi 

2003) . An additional six sources are reported 
with no detection in XMM-Newton references 
(IRAS 03256+3055; Preibisch & Zinneckcr 2002, 
HH211 MM; Preibisch & Zinnecker 2004, Ser- 
pens S68N, SMM2, SMM3, and SMM4; Preibisch 

2004) . Considering that these sources are deeply 
embedded, the lack of X-ray detection does not 
necessarily indicate that they do not have any 
X-ray-emitting activities, although they arc still 
difficult for X-ray as well as MIR investigations 
with current technologies. 
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